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The damping rate of hot giant dipole resonances (GDR) is investigated. Besides Landau damping 
we consider collisions and density fluctuations as contributions to the damping of GDR. Within 
the non-equilibrium Green's function method we derive a non-Markovian kinetic equation. The 
linearization of the latter one leads to complex dispersion relations. The complex solution provides 
the centroid energy and the damping width of giant resonances. The experimental damping widths 
are the full width half maximum (FWHM) and can be reproduced by the full width of the structure 
function. Within simple finite size scaling we give a relation between the minimal interaction strength 
which is required for a collective oscillation and the clustersize. We investigate the damping of giant 
dipole resonances within a Skyrme type of interaction. Different collision integrals are compared 
with each other in order to incorporate correlations. The inclusion of a conserving relaxation time 
approximation allows to find the T^-dependence of the damping rate with a temperature known from 
the Fermi-liquid theory. However, memory effects turn out to be essential for a proper treatment of 
the damping of collective modes. We derive a Landau like formula for the one-particle relaxation 
time similar to the damping of zero sound. 
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I. INTRODUCTION 

Giant resonances are high frequency, collective excita- 
tion modes of a nucleus. They can be identified as the 
collective motion in the nuclear volume and are found as 
a property of all nuclei. In particular in the last years 
the experimental and theoretical interest is focused to 
understand the width of such giant resonances [Q^^Oj. 
Although new experimental data are available at high ex- 
citation energy Jll| , p^ the theoretical description of the 
temperature dependence of damping rate is still a matter 
of discussion. 

The theoretical treatment of giant resonances can be 
roughly characterized by two approaches. The first one 
considers the finite nucleus and solves the random phase 
approximation ("RPA) equation by diagonalization of the 
Hamiltonian . In this approach the damping width 

can be extracted as an envelope of discrete excitation 
lines, sometimes called Landau fragmentation. However 
the introduction of the temperature remains difficult. 
The second class of approaches relies on the high den- 
sity of states and consequently uses a continuous model, 
mostly the Fermi hquid theory [|,0-|9[|l|-|l^ . Within this 
treatment, dispersion relations are derived whose solu- 
tions provide the energy and the width of collective ex- 
citations ||,|,|l3,0. 

The microscopic theory is mainly based on Vlasov ki- 
netic equations |]l8|-p0|| . The influence of correlations by 
particle-particle collisions is investigated using numerical 
solutions of Boltzmann-Uehling-Uhlenbeck (BUU)-type 
equations [^,^ . To get more physical insight into these 
simulation results, collective models based on the scaling 
theory are developed 123] . It turns out that the non- 



Markovian kinetic equation is necessary to get realistic 
values for giant monopole resonances These collec- 
tive models calculate the damping rate by an average 
procedure of the collision integral |24| . 

We follow another line. We will start from different 
kinetic equations and derive dispersion relations for the 
collective modes by linearization of the corresponding ki- 
netic equation. Instead of collective averaging, we solve 
these dispersion relations and obtain directly the influ- 
ence of correlations on the damping rate. We present in 
the paper different contributions to the damping of giant 
dipole resonances in a systematic way. 

Let us shortly outline our theoretical approach to col- 
lective resonances. The collective density oscillations are 
determined via the time dependence of the one-particle 
distribution function f{R,p,T) where the density reads 

n{R,T) = / f{R,p,T). The one-particle distribu- 
tion obeys the kinetic equation 



/(i?,p,r)-{i^ 



c/(i?,r) 

= /c[/(i?,p,T)], 



(1) 



where U is the mean-field potential and Ic corresponds 
to the collisional term. The Poisson brackets are ab- 
breviated as {a, b} = djiadpb — dpadnb. To get the 
linear response of the system to an external field Uext 
we linearize equation ffl) around a quasi-equilibrium 



fiR,P,T) = f{p) 
transformation T 



f Sf{R,p,T) and get after Fourier 
Lo and R ^ q 



i^5f{q,p,uj) 



■i—Sf{q,P,cd) 
m 



+ i{U'[n]5n{q,uj) + U,^)qdph{p) = h[5f{q,p,uj)]. (2) 
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Integrating over p the solution provides the polarization 
function Tl{q, to) 



niq,Lu) 



(3) 



with Sn{q,u!) = J Sf{q,p,uj). In the RPA approxi- 
mation which corresponds to the neglect of collisions we 
obtain 



nRPA('?, w) 



n"{q,Lu) 



1-Voll0{q,to) 



(4) 



with the standard form of polarization function 11° ( p^ ) 
and Vq = dnU[n\. The polarization function contains 
information about the collective excitation properties of 
nuclear matter. According to the denominator of (^) the 
relation to the dielectric function (DF) is given by 



e{q,Lo) ^1-Vo{no)n"{q,co). 



(5) 



For different collision integrals we will get different po- 
larization functions. The complex zeros ujq — — ij oi 
(H) determine the energy Q and width 7 of the collective 
excitation. ^From the DF one has the spectral function 
(structure function) via 



Ini- 



lme{q,uj) 



^) [Ree{q, iu)] + [lme(g, cj)] ' 



(6) 



which is important because it's structure reflects the col- 
lective excitations. Sum rules, e.g. |23] 



1 /■ , , 1 n^q' 
/ aw cjim— = 





(7) 



are an exact property of the spectral function. 

The paper is organized as follows. In Sec.|| we start 
with a generalized quantum kinetic equation which we 
rederived in appendix |^ using the Martin-Schwinger hi- 
erarchy for the real-time Green's function. To include 
memory effects into the kinetic equation we use the 
generalized Kadanoff-Baym Ansatz and obtain a non- 
Markovian kinetic equation. In Sec. II] we linearize these 
kinetic equations to get the DF for infinite h ot nu - 
clear matter in four different approximations: (Sec. [II A) 
collisionless Vlasov equation, (Sec. Ill B ) conserving re- 
laxation time approximation (Mermin approximation), 
(Sec p^II C ) dynamical relaxation time approximation (re- 



flecting memory effects) and (Sec. [[II D| ) the effect of 
density fluctuations on the potential. For these approx- 
imations we compare the damping rates and centroid 
energies of GDR which ar e the complex solution of the 
dispersion relation. In Sec. [HE we discuss these results 



together with the FWHM of the structure function and 
with experimental data. 



II. KINETIC EQUATION APPROACH 

Let us start from the kinetic equation in general form 
(A14) from appendix^ 



d_ 

df 



fw{p,R,T) 



dr 



dp' 
(2^ 



fwip',R,T) 



j:„{R+^,T)-^HiR-^,T) 



= Ic, (8) 



= dr 



J dr [{g>(p, R,T~ ^, r), E<(p, R,T ~ ^, -r)} 

- [G<{p,R,T~^,r),E>{p,R,T~^,-T)}] , 

(9) 

where the Wigner distribution function fw is connected 
to the correlation function fw{p, R, T) = {p, R,T,t = 
0). Further {, } is the anti-commutator of integrals over 
Wigner coordinates. 

Neglecting the collision integral Ic on the right hand side 
of (H) one obtains the collisionless quantum Vlasov equa- 
tion [26|. This leads to the Lindhard polarization func- 
tion (|l|). 

Next, we will consider binary collisions and will use for 
the self-energy in (O) the 1. Born approximation 



S<(p,i-t') 



d^p'd^p d^p' 
X W{pp',pp')[2TTf5{p + p'-p-p') 



(10) 



X G> (p', t' - t)G< {p, t-t')G<{p',t^t'), 

where W is the collision probability. We close equation 
(0) applying the generalized Kadanoff-Baym ansatz 



G<ip,R,T,T) 



-i / dte{p.,R.t) 



f^{p,R,T~^) (11) 



which gives a connection between the correlation func- 

> 

tions and the Wigner distribution. The quasi-particle 
energy e(p, R, T) in the quasi-particle picture is given by 
the solution of the dispersion relation 



P 

2m 



U{R,T) + ReJ:{p,e,R,T). 



(12) 



The resulting non-Markovian collision integral reads 
now 



d?P2<Pp'i 



I^niPuT) ^ J dr J '-^ 



W{piP2,P3P'i) 



[U+{T-T, T) + U- {T - T,T)] 
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X [/3(T - T)/4(r - t)(1 - /i(T - r))(l - /2(r - r)) 
- - r)/2(r - r)(l - /3(T - r))(l - ,U{T - r))] 

(13) 

where /,(r) = /(p„i?,T) (z = 1, 2, 3, 4), p4 = Pi +P2 -P3 
and the fuh time dependent propagator is 



C/± (T - r, T) = exp[±i / dt' A e{t')] 



(14) 



Ae are the time dependent quasi-particle energies 
Ae(i') = ei(t') + e2(i') - eM') ~ eiit'). If memory effects 
are neglected, equation (p^ becomes the usual Boltz- 
mann or BUU coUision integral Ib 



(27r)6 

X 5{e(j)i) + e{p2) - e(p3) - 4Pi) 

X [/3/4(l-/l)(l-/2)-/l/2(l 



/3)(l-/4)] . (15) 



The Bohzmann colhsion integral is modified in ( p^ ) 
by a broadening of the delta distribution function of the 
energy conservation and an additional retardation in the 
center-of-mass times of the distribution functions fi — 
f{pi, R,T — t). The first effect is connected with phase 
decay or spectral properties and responsible for global 
energy conservation |2^j2^. The second effect gives rise 
to genuine memory effects. The formation of correlations 
at such short time scales are discussed in lEgll. 



III. COLLECTIVE EXCITATION AND 
DIELECTRIC FUNCTION 

We now linearize the different derived kinetic equations 
and analyze their consequences to the damping of GDR. 
The standard RPA is first repeated in order to explain 
our analysis. 



A. Vlasov equation-collisionless Landau damping 

The linearization of the quantum Vlasov equation (|^) 
yields the RPA of the dielectric function which has the 
form of equation (||), but llo is the (complex) Lindhard 
polarization function pO| 



Ho (17, 7, g) = 2 



f{k)~f(k + q) 



(27r)3 Eu - Ek+q + (r! - ii) + iii 



(16) 



where E]^ = fc^/2m and rj is an infinitesimal small num- 
ber. Spin degeneracy has been accounted for. f2 and 
7 denote, respectively, the real part and the negative of 
the imaginary part of the frequency (a; = SI — 17). By 
consideration of a simplified Skyrme force [3l|] 



V = io(l + XQP„)5{ri - r^) + H8(rx - r2)S{r2 - rj) , 

(17) 

one obtains a mean-field potential for the neutrons Un 

US 



n{R, T) = to { (1 + y) {nn{R, T) + np(i?, T)) 
n„(i?,T) 



2^0 



+ ^ {(n„(i?,T) -|-n,(i?,r))2 - n2(i?,T)} (18) 

and Up is given by an interchange of Up and 7i„. /,From 
(^) and (H) we can read off the effective particle-hole po- 
tential Vq for isovector modes 



Un-Up = 2Vo [Sun - (^rip] 



Vo(no) 



to 



(jln + Up) . 



(19) 



The parameter ioi ^3 a-nd xq were fitted to reproduce 
the binding energy (Ej A = —16 MeV) at the saturation 
density (jiq — 0.17 fm""^) of nuclear matter. For the GDR 
the wave vector g is estimated according to the formula 



9 = 



2R 



(20) 



In this model neutrons and protons oscillate out of phase 
inside a sphere of the radius R (=nuclear radius). Fo- 
cusing our interest to the nucleus ^ospj-, (i20gjj-j .^^g 
q K. 0.23 fm^i (9 « 0.277 fm^i) with R = 6.7 fm 
{r = 5.6 fm). 

The dispersion of collective excitation is now computed 
from the zeros of the complex dielectric function (DF), 
Eq. d). 



Re e(ri — 27, g) + Im e(Sl — 27, g) = 



(21) 



where 7 gives the Landau damping of the collective exci- 
tation. An approximate solution of this RPA dispersion 
relation is possible if the damping Ime is small . Then 
one can linearize the collective excitation spectrum 

Ree(g, f7) + (cj - f7 + i7)a^Ree(g, ^) + ilme(g, 17) = 

(22) 



which leads to 



Ime((7, ri) 
5nRee(g, 17) 



(23) 



and r2 the solution of Ree(g, SI) = 0. This is, however 
only justified for small values of the damping 7. The 
correct procedure is to carry out the analytical continu- 
ation of the DF into the lower energy plane. Performing 
the integration one can express the DF (H) with ( [l^ ) by 
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the dimensionless variables x = q/\/2mT, z 
and ^ = — 7/2T in the form 



F{z) 



F{z) = hi 



1 + i^cxp- (z - I)" 

1 + (TCXp - + f )" 



where (24) 
(25) 



c = ypXTrr? jix^ and = e^/-^ denotes the fugacity. Fol- 
lowing Landau's contour integration |p5| the result of 
the analytic continuation of the DF ( p4[ ) with the pole 

zo = 2* — ?^ is 



e-[x,z\i) 



Voc 
2x 



+ CXD 



dz 



—00 



dz- 



F{z) 



F{z) 



inFiz*) C = 

+ i2T:F{z* -i^) C > 



dz- 



F{z) 



(26) 



The explicit expressions of the real and imaginary parts 
of the analytical continuation of the DF are given in ap- 
pendix p. For complex values of z there are also poles of 
the function F{z) in ( p5[ ) which require a separate inves- 
tigation. They are located at 



(27) 



where dn are the discrete fermionic Matsubara frequen- 
cies 




i9„ = /i- i(2r7. + l)7rr; 7i = 0,±l,±2, . 



(28) 



If these poles do not agree with the poles zq of the denom- 
inator of equation (M) there is no contribution to the in- 
tegration due to the fact that Kes[F{z)/{z — zq); z„] = 0. 
Because of the infinite residue the remaining case zo = z„ 
is found to be singular and will be discussed elsewhere. 

Fig. |l| shows the zeros of the real and the imaginary 
part of the DF equation (|2^) in the complex plane at 
zero temperature (left hand side). The right hand side 
represents their position in the pair excitation spectrum 
which is bounded by the lines w± = q/2m(q± 2pp). For 
increasing wave vectors q we have found three special 
cases: 
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FIG. 1. Zeros of the analytic continuation of the complex 
DF (Lindhard) for different wave vectors at T = 0. LHS: 
collective excitations (GDR) in nuclear matter correspond to 
the crossing of the lines Re e = and Im e = (marked by 
thick dots). RHS: pair continuum Im e 7^ and undamped 
(marked) region Im e = 0, 7 = with uj± — q/2m(q ± 2pf)- 

Top (q=0.23 fm~^): There are two collective excita- 
tions which correspond to the crossing points of the zeros 
of Re(e^) and Im(e'^). Their resolution in the pair con- 
tinuum (RHS) yields an undamped collective excitation 
LU02 which lies outside the pair continuum (marked re- 
gion, Im(e) = 0, 7 = 0). The collective excitation wqi 
lies inside the pair continuum and is therefore damped 
(Im(e) 7^ 0, 7 7^ 0). From the DF one can now calculate 
the spectral function (^). In Fig. || one recognizes a S- 
shaped peak at the frequency of the undamped excitation 
LU02 and a broader peak corresponding to the damped col- 
lective excitation ujqi. 

Middle (q=0.338 fm^^): There arises only one critical 
collective excitation cocrit for a critical wave vector qcrit- 
The collective excitation ujcrit hes inside the pair contin- 
uum {tOcrit ~ i^+) and is damped. 

Bottom (q=0.4 fm~^): Since there is no crossing of the 
zeros of the real part and the imaginary part of the DF, 
no collective excitation can occur. 
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FIG. 2. Real part (solid line), imaginary part (dashed line) 
and spectral function (dotted line) of the Lindhard DF for the 
wave vector q — 0.23 fm~^ at T = corresponding to Fig. |l| 
(top). 

In Fig. ^ we summarize the above results considering 
the entire dispersion of the collective excitation at T = 0. 
We see that nuclear matter at zero temperature has a 
region where two collective modes are excited (Fig. ^ 
(A)). The mode r2o2 goes outside the pair continuum and 
is undamped, whereas the mode f2oi propagates inside as 
a damped one. 



Beyond the coincidence of both modes into the critical 
point ^critiQcrit) ~ i-^+iQcrit) there are no further collec- 
tive excitations. The related damping rates are shown in 
Fig. H (B). We observe that above a critical wave vector 
we cannot find collective excitations which represents a 
pure quantum effect. This critical wave vector Qcrit is 
determined by the used interaction Vq. Now we can link 
the minimal interaction required for collective oscillation 
with the mass number by equation ( pO| ) . The relation be- 
tween minimal interaction Vq^^^ and mass number Acrit 
(wave vector qcrit) is plotted in Fig. ^. As a result we 
find the following fit 

Vo^„, {Ac„t) = (a + MeV fm^ , (29) 

where a = 151.9 and b — 2997.4. The marked area Fig. 
^ designates a region where no collective excitations ex- 
ist. One sees that a certain interaction strength is neces- 
sary in order to have collective oscillations. The physical 
meaning is that for light nuclei with given interaction 
strength the sound velocity couples to the single particle 
motion more strongly. This fact is reflected in the criti- 
cal wave vector, where the collective mode enters the pair 
continuum (cp. Fig. 
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FIG. 3. (A) Dispersion of the collective excitations at 
T = of Fig. hi and their position in the pair continuum. The 
mode no2(<?) (dotted line) is undamped and Qoi{q) (dashed 
line) is damped. The insert shows an enlarged view of the 
pair continuum where the marked region corresponds to the 
undamped region (lm=0). (B) Corresponding damping rates. 



temperature T=0 MeV 
density n„=0.17 fm"° 




100.0 150.0 

mass number A„ 



FIG. 4. According to Fig. |l| (middle) calculated critical 
potentials Vo^^^ for critical wave vectors qcrit which are con- 
nected with the mass number A via (|2o|). The marked area 
corresponds to a region where no collective excitations exist. 



Since this model relies on the simple scaling law (|2C 
the equation (|2^) can only be a qualitative estimation. 
For finite nuclei we expect different parameters a and b. 

Considering now nuclear matter for finite temperature 
we present in Fig. ^ the zeros of the real and imaginary 
part of the DF compared with the T = result. We find 
with growing temperatures a continuous deformation of 
the lines Re(e^) — and Im(e'^) = 0. The crossing points 
(black marked) correspond again to collective excitations 
and are always damped. We have additional zeros (light 
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marked) in contrast to the zero temperature result. They 
correspond to the complex poles of equation (|2^) and 
can be identified as single particle excitations. Taking 
the entire dispersion relation for different temperatures 
(Fig. ^) we observe higher damping rates of the collec- 
tive excitations as the temperature increases. This effect 
of temperature is also visible in the calculated spectral 
function (|) of Fig. 0. 



Growing temperatures lead to a broadening of the peak 
and reaching a temperature of 6 MeV (solid line Fig. 
^ only small collective effects remain observable. The 
corresponding centroid energy is increasing with higher 
temperatures. The further inclusion of correlations will 
lead to a decreasing of the centroid energy for higher 
temperatures, which has been demonstrated in p7|,|31|. 
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FIG. 5. Collective excitations (thick dots) corresponding 
to the zeros of the analytic continuation of the complex 
DF (Lindhard) at different temperatures compared with the 
T = result (top or Fig. |l| (top)). 
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FIG. 7. Imaginary part of the inverse Lindhard DF (spec- 
tral function) for different temperatures and the wave vector 
q — 0.23 fm~^. (This figure corresponds to the results of pl|). 

In Fig. ||we see now for the two nuclei ^^''Sn and ^ospb 
(dotted hues) the calculated (Landau) damping rates of 
GDR over the nuclear temperature. These damping rates 
are compared and discussed with other approximation in 
the next sections. 



B. BUU equation-relaxation time approximation 
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FIG. 6. Collective excitations at T = 2 MeV (sohd line) 
and T = 6 MeV (dot-dashed line) of Fig. |^. (A) Dispersion 
of the collective excitations, (B) damping rates. 



In order to consider collisions as an additional damping 
effect of GDR we start first from (^ with the Markovian 
Boltzmann (BUU) collision integral ([l^ ) 1b- 

Within the standard treatment ||35|| we linearize the 
collision integral around a homogeneous equilibrium dis- 
tribution 

SMp,R,T)^f,{p,R,T)^f^{ep) 

= Mp,R,T)^f^ (30) 

where should be first the global Fermi distribution. 
Later on we will use a local equilibrium distribution in 
order to ensure conservation laws. The linearized colli- 
sion integral reads 

X [/lV2°(l-/3)(l-/4)], (31) 



6 



where Ae = £1+62 — £3 — £4 and A<I> = $i + $2 — 3>3 — $4. 
Neglecting the backscattering terms ^2 — ^3 — $4 we 
obtain the relaxation time approximation 



with 



f-f _ Sf 



(32) 



tb{pi) 



(27r)6 

/2°(l-/3°)(l-/4") 

1-/1° 



M^(12,34)^(Ae) 



(33) 



Furthermore we will use a thermal averaged quantity. 
This procedure is introduced in appendix ^ and we get 
from das'' 



1 

TB 



1 

no 



(Ppi(fip2d^P3 



(27r)9 

X [/i°/2°(l-/3°)(l-/4°)] 



W {piP2,P3Pa) 5{/\e) 



(34) 



with nu = mppT/{2TT^). For temperatures which are 
small compared to the Fermi energy we follow the well 
known methods of the Fermi liquid theory |^,^ and get 
from (11) 



1 

TB 



2TTm% / W{d,cl>) \ 
riB \ cos6/2 I " 



B 



where 



(...) 



27T 



2tt 



(27r)s 



de sin 



cos 61/2 



(35) 



(36) 



The calculation of the Fermi integral Ig is performed in 
appendix ^ with the result 



27r2 
3 



.y3 



(37) 



Using the collision probability in Born approximation we 



introduce the spin-isospin averaged cross section ^ 



1^(61,0) = 27r 



m 



da 

dn 



with the spin-isospin degeneracy g — A. 
mean cross section ^ ~ ^ we get finally 



1 

TB 



-T ma , 



(38) 

Assuming a 
(39) 



where a is the nuclear cross section (for numerical calcu- 
lation we use a — AO mb as a spin-isospin averaged cross 



section of ct„„ , app and cr„p 




1.0 2.0 

temperature T in [MeV] 



FIG. 8. Theoretical damping rates of GDR for ^^"Sn (grey) 
and ^'^^Pb (black) as a function of the nuclear temperature T. 
They are the complex solution of the dispersion relation for 
different approximations: collisionless (Lindhard) DF (dot- 
ted line) , Mermin DF (Boltzmann) (dot-dashed line) , Mermin 
DF (memory) (solid line) and Mermin DF (density fluctua- 
tion-!- Coulomb interaction) (dashed line). 

Based on this result one can use the relaxation time 
approximation to include collision effects into the DF. 
This leads to a pure replacement lo ^ lu + i/T in the col- 
lisionless Lindhard DF (^6|) and is known to violate the 
sum rule (0). Therefore we apply an extended approach 
using a conserving relaxation time approximation follow- 
ing Mermin pof . The local equilibrium distribution /"^ is 
characterized by 



cxpL 



li-Sn{R,t))] + l 



(40) 



where the local chemical potential Sfi is determined by 
the density conservation resulting in 6p,{q,LLj) = j^r^^- 
The linearized kinetic equation M) reads then 



q-p 



Sf{p, q,w) 



[f{p-q/2)-fip + q/2)] SV{q,u;) 



1 



Sf{p, q,(^) 



df{p) 5n{q,uj] 



TB \ dep no((?, 0) 

and the resulting (Mermin) DF is 



(41) 



{q,u) = 1 + 



1 



UJTb 



[e{q,uj + i/TB) - 1] 



UJTB 



[e{q,uj + i/TB) - 1] 



(42) 



Here e is the collisionless Lindhard DF (to calculate the 
zeros of e*^ one has to use the analytical continuation of 
Lindhard DF equation (p6|)) and tb the relaxation time 
(^) . Further conservation laws can be incorporated Eol . 
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In Fig. ^ we compare the spectral functions com- 
puted with the Lindhard DF (dashed Une) and Mermin 
DF (sohd hne), respectively. We observe an additional 
broadening beyond Landau damping in the width of the 
Mermin DF due to collision damping in this approxima- 
tion. In comparison with Fig. Q we see further that the 
centroid energy of Mermin result (solid line) is shifted to 
lower values with increasing temperature, while the cen- 
troid energy of the Lindhard result (dashed line) is shifted 
to higher values. The result for the damping rates 7 of 
the GDR of i^OSn and ^°^Ph is shown in Fig. | (dot- 
dashed lines) which are determined by the zeros of the 
Mermin DF e*^(ri, 7) = 0. The calculated damping rates 
in this approximation present an improvement compared 
to the collisionless results (dotted lines) . We observe the 
typical behavior of the damping known from experi- 
mental data (Fig. |l^) but we do not get any contribution 
at zero temperature. Now we shall employ a better ap- 
proximation of collisions in the next section. 
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FIG. 9. Spectral function (|) with the Lindhard DF 
(dashed line) compared with the Mermin DF for different re- 
laxation times: Boltzmann relaxation time tb (solid line) , dy- 
namical (memory) dependent relaxation time Tm (dot-dashed 
line). 



C. Non-Markovian equation-dynamical relaxation 
time 

Next we investigate the kinetic equation with the col- 
lision integral (tt3[) taking into account memory effects. 
The linearization of the collision term ( |l3|) can be 
performed in the same way as above with an additional 
linearization with respect to . 
gives 



Expanding U in (llJ 



T-r 



SU^ (T - T, T) = {T)[l±i J dt' A Se] , (43) 

T 

where 

A6e = Aeit') - a£^ and U^ir) = exp(±z(A|lr)). 



The linearization of the colHsional term ( [13| ) reads now 

U+{r) + U^{r){SF,^5F2) 
f dU+{T -T,T)- 6U- (T - T, T) 



(44) 



where 



- inm - - /2) - /?/2"(i - - /4°)]- 

(45) 

The calculation of 5Fi with the help of (pO|) results into 



[fl^l + /2$2 - (1 - /3°)'i>3 - (1 - /4°)$4] 



T 



(46) 



and 5F2 is obtained from equation ( [iq) by interchanging 
/ ^ (1 — /). According to equation (^ we perform a 
Fourier transformation T — > a; of the linearized collision 
integral (E3) with the result 



/™(pi, w) 



1^(12, 34) 



D+{lo) {SFi -SF2) + A 6e 



OJ 



(47) 



where D±{iL;) = 6{lj + Ae) ± 6{lj - AeV For 6e in (|4 
we obtain from the dispersion relation ( [l2|) 



U\no)5n{R,T) 

(48) 



where 5n{R,T) — '^p5f{p,R,T). In the following we 
neglect self-energy effects (S = 0) and restrict our consid- 
eration to quasi-particle energies in mean-field approxi- 
mation Se{pi^R,T) = U'ino)Sn{R,T) with U' = Vo from 
equation Therefore we see that A6e — 0. This 

would be not the case for momentum dependent poten- 
tials. 

One can now introduce a thermal averaged relaxation 
time r,„ analogous to (p4) 



1 



1 



1^(12,34) 



X |z?+(c.) (/0/3°/4° - f^flf!) } , (49) 

where / = (1 — /) and tt-d is specified in appendix JQ. 
Comparing the Markovian result equation (|3J) with ( |49[ ) 
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we observe a replacement of the energy conserving factor 
S{e) with S{uj + Ae) and S{uj — Ae). The energy huj of 
the collective excitations are now included into the energy 
conservation. We can interpret this effect as a coupling of 
the collective modes to the binary collisions. The collec- 
tive boson is absorbed or emitted when two particles are 
colliding ]3^ . Restricting the relaxation time t™ again 
to temperatures which are small compared to the Fermi 
energy we follow (Isl) and get from (p9|) 



1 



2Wp/ /W{d,(l>)\ 



/ y \ jf / X 



(50) 



where the Fermi integral ll^ is done in appendix ^ with 
the result 



rj-io 



— vr 
3 



2 vr 



The analytical expression for finally reads 



1 



T„i(w) 



1 

TB 



4 \7tTJ 



(51) 



(52) 



where tb is the relaxation time (p9| ) which corresponds 
to the Markovian limit. The expression (^2|) for the re- 
laxation time Tm contains an additional (dynamical) con- 
tribution (~ w^). It guarantees that also at zero temper- 
atures the collective mode can be considered as a self- 
propagating one which has a finite damping. We obtain 
a relaxation rate which is the Landau damping rate of 
zero sound except the factor of 3 in front of the fre- 
quencies. 



The form (|52| ) is similar to the result Ayik derived in 
p8[ for the collective damping rate. In contrast to this 
result we have here the single-particle relaxation time. 
As long as we have only momentum independent mean 
fields, the linearization of the mean field propagator does 
not reduce this to the standard Landau result as claimed 
in the errata of |^^, instead the factor of 3 counts here 
(see comment after equation (p8|)). 

One may argue that this factor 3 comes from the ther- 
mal averaging we employed here. We can repeat all cal- 
culations without thermal averaging (cp. appendix Dl) 
with the momentum in the relaxation time at the Fermi 
level. The resulting dynamical relaxation time reads then 



1 



3 1 

4ts 



(53) 



Compared with (^2|) one sees no factor 3 appears in front 
of the frequency, however a factor 3/4 in front of the re- 
laxation time 1/tb- Therefore we obtain for T = the 
same frequency dependence in both cases. In other words 
the zero sound damping is not affected by averaging. The 
temperature dependence is slightly flatter without ther- 
mal averaging. 

In the Mermin approximation e^^ of (^) we have to 
replace the static relaxation time tb with the dynamical 



one rm(w). Then we can again compute the damping 
rates 7 for the GDR of ^^°Sn and ^°*Pb by searching 
for the zeros of the Mermin DF e^'^{fl,j) = 0. The 
results are shown in Fig. || (solid lines). We observe 
a T^— dependence of the damping which is due to the 
low temperature expansion of the relaxation time t,„, 
but the temperature increase is too flat compared with 
the Boltzmann result (dot-dashed line) . The inclusion of 
memory effects now induces a finite width at low temper- 
atures. This reflects the fact that memory accounts for 
zero sound damping. Next we investigate the influence 
of density oscillations on the potential itself. 



D. Contribution of density fluctuations 

As shown in ^ the contribution of density fluctuations 
can have an remarkable effect on the damping rate. The 
derivation of the kinetic equation including the density 
oscillations leads to Lenard-Balescu type of collision inte- 
grals which have the same structure as (|l^) but with 
a dynamical potential 



(54) 



replacing the static one. Here the DF renormalizes the 
potential. Then we can proceed as described above and 
linearize the collision integral, etc. This would lead to a 
rather involved integration. We simplify the treatment 
by assuming that the transfered momentum during a sin- 
gle collision is small. Since for small q we have for the 
Lindhard DF 



\une{q,Uj) = l-'^+Oiq^) 



(55) 



with the free sound velocity — Von/m, we get from 
(1) 



lim e^\q,uj) = 1- 



2 2 
c q 



q^O 



lo{lo -\- i/r) 



(56) 



Inspecting ( ^ ) we see that the frequency argument of 
( p^ ) has to be taken at ei — e'l = 9^^^Tw « if we fix 
the final state of scattering as the state which matters for 
kinetic processes and use vanishing transfer momentum 
q. Further we observe that \e{q,u})\ = |e(q, — w)| such 
that we have a common factor in ( |49| ) in front of 
of \/\e{q 0,uj)\'^. The further procedure is as above 
described. We expand for low temperatures and observe 
that for the angular averaging holds 



q 



cos 61/2 



2pF sin(6'/2)sin(0/2) 
cos6'/2 



(57) 



Using ( pq ) the resulting prefactor which renormalizes 
(M) reads 
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dx \e ' {2ppx, uj)\ 



dx- 



with 







(58) 



(59) 



w(^ + ^) 

and the complex frequency uj ~ Q ~ i"f. With tabulated 
integrals |^ we get from (|58| ) finally an analytical ex- 
pression 

1 Ik 



tlb{^,i) T,„(ri, 7) 2 sin 2a 

sin a arctanh(2fc cos a + fc^) 



cos a arctan | 



1 



2fcsinQ! 



(60) 



where a = ^ arccos (Rez/jzp) and k — 1 / \/\z\. For 
our chosen situation and potential the prefactor m^) in- 
creases the damping rates only by 0.2 MeV. However, 
we see that the renormalization of the potential by den- 
sity fluctuations which are in turn determined by the lin- 
earization of the kinetic equation can lead in principal to 
a remarkable change in the dispersion relation if we come 
close to the instability line |^ . Our model potential ( pj| ) 
has no instability in the parameter range. The proce- 
dure here means that density fluctuations are caused by 
interactions or correlations, however these density fluc- 
tuations renormalize the potential, i.e. their cause. So 
we have a complicated feedback of correlations to the 
fluctuations and so forth. 



1. Inclusion of Coulomb effects 

We expect a more pronounced effect if long range in- 
teractions mediate collective oscillations. The Coulomb 
interaction leads to an additional contribution for the 
proton meanfield (^8|) 

Ane^ , , , , 

—np{q) , (61) 



so that the resulting dispersion relation will be a matrix 
equation As above, the damping rates 7 are again 
the complex solutions of this dispersion relation and plot- 
ted in Fig. H (dashed lines). For small temperatures we 
find an increasing of the damping rates compared with 
the rates of the memory dependent collision approxima- 
tion ( [5^ ) (solid lines). This comes from the fact that we 
have density fiuctuations at zero temperature caused by 
the Coulomb interaction. For temperatures larger than 
3 MeV the binary collisions dominate and the behavior 
follows the memory collision approximation. 



E. Comparison with the Experiment 

It is interesting to point out an advantage of the Mer- 
min polarization function. Therefore we compute the 
power spectrum of the mode, i.e. the energy rate per 
time which is expended on the motion of the collective 
mode. This power spectrum is connected to the structure 
factor S{q,ui) = J2 \<f \Vo\0> \'^S{huj - Ef + Eq) by 
/ 



27r 

P{q,i^) = —S{q,uj)huj 



(62) 



which is just Fermi's Golden Rule. The structure factor 
itself is given by the dielectric function S = Vo/7rIme^^. 
Using (pq) we arrive at an expression for the power spec- 
trum of 



P{q,uj) = 2Vbw Im 



(63) 



with ll>q = Vgnq^ /m. The second line is derived with the 
assumption that t is real, which is not true for dynamical 
relaxation times (memory effects). We rederive by this 
way just the classical Lorentz formula which describes 
the energy rate expended on the motion of a damped 
harmonic oscillator driven by the external force x + x/t + 
luqX — 2Voujq cos{ujot) averaged over time. 




temperature T [MeV] 

FIG. 10. Temperature dependence of the centroid energy 
Egdr and of FWHM for GDR mode in ^ospb (grey lines) 
compared with the solution of the complex dispersion relation 
(pl[), Retj = Q and Imuj — 7_(dark lines), using the memory 
dependent relaxation time ( |52[ ) in the Mermin DF (^). 

Naturally that ( |63| ) leads to a Breit Wigner form near 
the resonance energy loq with the full width of half max- 
imum (FWHM) of r = l/r 
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P{q,Lo) = VoiU^ 



(c.-c.o)2 + (E)2- 



(64) 



The damping rate in classical approximation ( |56| ) (long 
wavelength) is given as 7 = 1/(2t). We recognize that 
the FWHM is just twice the damping rate F = 27. This 
has been recently emphasized It has to be stressed 
that the experimental data are accessible by FWHM. To 
extract the FWHM from the structure function (||) we 
used the Mermin approximation e*^ (^) with the dy- 
namical (memory) dependent relaxation time r,„ of ( |52[ ) . 

In Fig. |l^ we have plotted the temperature dependence 
of FWHM (grey-dashed line) of the structure function 
(cp. Fig. ^ and the corresponding centroid energy 
(grey-solid line) of the GDR mode in ^°*Pb together with 
real part (dark-solid line) and imaginary part 7 (dark- 
dashed line) of the complex solution of the dispersion 
relation, respectively. We find an approximate relation 
FWHM « 27 (grey-dashed and thin dark-dashed lines) 
which holds in the temperature limit, where FWHM < 
Eg DR., as has been discussed in 

We emphasize that the centroid energy Egou of the 
structure function agrees with the real part 17 of the com- 
plex solution of dispersion relation up to higher temper- 
atures (thick-dark and thick-grey lines). The small de- 
viation from the exact Breit Wigner form ( |6^ ) is caused 
by the memory effects resulting in a frequency dependent 
dynamical relaxation time t„i which is now complex. 




temperature T in [MeV] 



FIG. 11. Experimental damping rates of GDR for ^^°Sn 
and Pb (i20sn from and ^ospb from [||) compared 
with the FWHM of the structure function as a function of 
the nuclear temperature T. 

The FWHM is the observable which allows us to com- 
pare with the experimental data. In Fig. O we have 
plotted the FWHM of the GDR modes in ^Sn (grey 
line) and ^°^Pb (dark line) as a function of temperature 
compared with experimental data. They coincide within 
theoretical limits if we consider that the experimental 
values are fits of a Breit Wigner form itself. 



In Fig. 12 we compare different approximations calcu- 
lated in this paper with experimental data versus mass 
number for T = 0. In (A) we have plotted the cen- 
troid energy Vl — Egdr- We observe that the suc- 
cessive inclusion of collisions with the memory depen- 
dent relaxation time t„i (^) (thin solid line), density 
fluctuations with the relaxation time tlb ( |60| ) (dashed 
line) as well as the inclusion of the Coulomb interaction 
(dotted and dot-dashed line) reproduces the experimen- 
tal values (diamonds) increasingly for the mass numbers 
80 < A < 210. The inclusion of only Coulomb effects 
(dotted line) slightly overestimates the data. This overes- 
timate is compensated if we add the density fluctuations 
(dot-dashed line). The corresponding FWHM, which are 
shown to be roughly twice the imaginary parts 7 of the 
complex dispersion relations are plotted in Fig. [l^(B). 
The inclusion of collisions brings the curve towards the 
experimental values. The improvement by inclusion of 
Coulomb effects and density fluctuations is small. 

temperature T= MeV 
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FIG. 12. The experimental centroid energy Egdr (A) and 
the experimental data of the damping widths 7 (B) for GDR 
psl vs mass number together with different approximations. 



IV. SUMMARY 

We have presented a systematic microscopic descrip- 
tion of collective excitations in hot nuclear matter by a 
kinetic approach and computed damping rates of giant 
dipole resonances (GDR). We have rederived a general- 
ized quantum kinetic equation employing the real-time 
Green functions technique. Within Hartree approxima- 
tion one derives the coUisionless Vlasov equation which 
in the linearized limit leads to the known Lindhard DF 
in RPA. Applying the complex integration method pro- 
posed by Landau we calculated the damping width (Lan- 
dau damping) of GDR in different nuclei at finite tem- 
peratures. 

The experimental data are calculated via the FWHM. 
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Landau damping rates underestimate the experimental 
data of ^^°Sn und ^°*Pb. Considering collisions in first 
Born approximation the (Markovian) kinetic equation is 
of Boltzmann or BUU type. Using a relaxation time ap- 
proximation we have incorporated the collision effects 
into the DF which leads to the Mcrmin DF. The cal- 
culated damping rates of the GDR show a typical T^- 
behavior but do not reproduce the experimental values 
yet. Only memory effects which are included into the 
collision term of the kinetic equation via retardation of 
the distribution function improve the theoretical damp- 
ing rates. This non-Markovian relaxation time approxi- 
mation leads to a dynamical relaxation time which differs 
from the Landau result of the damping of zero sound. 
Within the low temperature limit the theoretical damp- 
ing rates of the GDR in ^^"Sn und ^°^Pb are improved 
but the temperature increase is flatter than the exper- 
imental findings. A renormalization of the potential by 
density fluctuations increase the damping rates by « 7%. 
These density fluctuations together with Coulomb ef- 
fects reproduce the centroid energy for the mass numbers 
80 < A < 210. If we neglect Coulomb interaction we un- 
derestimate the centroid energy for heavier elements. 

We consider the solution of complex dispersion rela- 
tions as the physical mode of collective excitation in the 
system. In contrast, the experimentally observed collec- 
tive excitations are extracted from FWHM. The calcula- 
tion of the width (FWHM) of the GDR modes, which is 
proportional to imaginary part of the inverse DF, leads to 
results which are comparable to the experiment. Mem- 
ory effects change the shape of these structure functions 
compared to the usual Breit Wigner form. Since exper- 
imental data are fitted with the latter one, we cannot 
expect to reproduce the experimental data completely. 
Especially the temperature dependence of 7 as solution 
of the dispersion relation shows too small increase with 
temperature. This may be a hint that shape fluctuations 
become important |Q. In it was stated that 7 in- 
creases linearly with temperature, caused by the coupling 
of surface modes. The incorporation of shape fluctua- 
tions will be considered in a forthcoming paper. 
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Green's functions introduced by Kadanoff and Baym [g6| . 
Considering a system of Fermions which interact via the 
potential V we start with the Hamiltonian 



H 



2 m 



+ i y" d\id?r2^\^lVi2{ri ~ r2)^i^2, (Al) 

where the different numbers, i.e. 1 = (r, i, s,...), denote 
the one-particle variables (space, time, spin,...). The an- 
nihilation and creation operators for Fermions obey the 
commutator relations 



[*l,*2]+=0, [^-1,*^]+ =(5i2. 



(A2) 



The correlation functions are defined by different prod- 
ucts of creation and annihilation operators in the Heisen- 
berg picture 



G>(l,2) = (*i*|), G<(l,2) = (vl/tv|/i) 



(A3) 



Here the (...) denotes the average value with the un- 
known quantum-statistical density operator p. The 
causal Green function is given with the Heaviside step 
function O by 

G(l,2) = e(ti -t2)G>(l,2)-e(t2-ti)G<(l,2) (A4) 

It is furthermore useful to introduce the retarded and 
advanced Green functions according to 

G^(l, 2) - - iQih - t2)[G> (1, 2) + G< (1, 2)] 
G-4(l,2)= ze(i2-ii)[G>(l,2) + G<(l,2)]. (A5) 

Using now the equation of motion of the creation and 
annihilation operators we can derive kinetic equations. 
Applying the equation of motion for the field operators 
in the Hcisenberg picture one gets an equation of motion 
for the causal Green functions fE6l 



zA + Zi)Gr(l,l') = 
oti 2m / 

(5(1 - 1') + j d2 V{1 - 2)G2(12, 1'2+) . (A6) 
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APPENDIX A: GENERALIZED KINETIC 
EQUATION 

We shortly sketch the derivation of the general ki- 
netic equation (H). We use the dialect of the general- 
ized Kadanoff and Baym formalism developed by Lan- 
greth and Wilkins for the nonequilibrium (real-time) 



In this so-called Martin-Schwinger hierarchy the one- 
particle Green's function couples to the two-particle one, 
the two-particle Green's function couples to the three- 
particle one, etc. A formally closed equation for the one- 
particle equation can be reached with the introduction of 
the self-energy 

y"diE(l,T)Gi(l,l') = yd2U(l-2)G2(12,l'2+) (A7) 



where the integration contour c turns out to be chosen 
as Keldysh contour in order to meet the requirement of 
weakening of initial correlations 
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lim G2(12, 1'2') 

{tl—t2—t[—t'^)—>- — OG 

= G(l, l')G(2, 2') ± G(l, 2')G(2, l') 



(A8) 



which means that asymptoticaUy the higher order corre- 
lation s sh ould be f actorized for infinite past. One gets 
from (A6) and (A7) the Dyson equation 



+ / dl (E(l,l)Gi(l.l')-S'=(l.l)Gf(l,l')). (A9) 



Writing (A9) into the more compact operator notation 
gives 



Gq ^Gi = 1 + SGi 



(AlO) 



where Gg ^ is the inverse of the Hartree-Fock Green func- 
tion 

G„-i(l, 1') = (^z^ + Yl^ 3(1 1') - Ehf(1, 1') . (All) 



We apply now t he La ngreth-Wilkins rules [Q onto the 
Dyson equation (AlO), which describe the way to get the 
correlation functions (A4) and retarded functions (A5) 
from causal ones. We get the equation of motion 



(Gn ' - S"") G< = S<G^ 



(A12) 



Subtracting the adjoined equation the general kinetic 
equation reads l47| 



I (G^I]< - I]<G^) - i (I]^G< - G<S'4) , (A13) 

which was derived first by Kadanoff and Baym P6[. For 
the time diagonal case t\ = ty = t and using ([A5[ ) in 
( A15| ) we finally get 



-t[Go\nt,n,t),G<{nt,nd)] = / d¥i / dt 



G>{nt, 1)S< (1, ryt) + S<(rii, 1)G> (T, ryt) 
-G<{nt, T)S>(1, ryt) ~ E>(rit, 1)G<(1, ryt)^ , (A14) 
which is the kinetic equation (pf) and (H), respectively. 



APPENDIX B: EXPRESSION OF THE 
DIELECTRIC FUNCTION 

In this appendix we give the explicit expressions of the 
real and imaginary parts of the analytical continuation 
of the complex dielectric function e'^(g, n, 7) in Eq. (p6|) 



Ree^(<7,f],7)- 



Ime^(g,r!,7) 



eRe{q,^) , 7 = 
ReeRe{q,fl- 17) 

-2Ime/„(g, fi- 17) , 7>0 

ReeReiq,^- il) , 7 < 



e/m(q, rj) , 7 = 

+2Ree/™(g,fJ - 17) , 7 > 



(Bl) 



(B2) 



where we redefined e"^ in terms of the retarded Lindhard 
DF e{q,uj) — eRe + i^im (P^. For nonzero temperature 
we find 



e_Re = 1 



Vom 
TT^q 



dkkfpik) 0(y+)-(/.(y_) , (B3) 



VoTm^ 



In 



nq 



l + exp^('^-(^-)'/2m) 



l + cxp'^f^-t^^^'/^™) 



(B4) 



with y± = B%/k, B% = Lom/q ± q/2, and ^(y) = 

ln|(l-y)/(l + 2;)|. 

Considering now complex frequency u — fl — i^ we find 
for e_Re(g, 0,7) = Ree_Re -|-Ime_Re 



Re €Re = 1 - 



X < In 



27r2g 



dkkfpik) 



(B5) 



(fc2 _ 



(2/ck)2 



-In 



[(fc+(B^^)2-K2]" 
(fc2 - {B'lf - K^f + {2kKf ' 
[{k + B^}y - K^]^ 



(B6) 



and 



^ / dfcfc/F(fc) 



(B7) 





arctan 



2kK 



arctan 



(fc2 - (B^^)2 - k2) 
2fcK 



(fc2 - (BL^)2 - ^2) 



(B8) 



where k = ^m/q. 

The expressions of Re e/m and Im e/™ are 



Re e/„i = 



2Trq 



In 



and 



with 



Im e/m. — arctan — 

■nq \ac 



(B9) 



(BIO) 
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(7 = 1 + 2a_ cos{B^Ki) + at 

(Tc = 1 + a+ cos{B^2ki) + a_ coa{B^Ki) + a+a_ cos{—Kiq) 
(7s = 0,+ sin(i?^Ki) — a_ sin(i3^Ki) + a+a_ sin(— Kiq). 

(Bll) 

The coefficient a± is defined by 



where /° = (e^ + = 1 - /° and Ax = + - 

X3 — X4. The integral (Dl) can be done exactly using the 
standard identity over Fermi functions^ [3^ ] and we find 



2^3 



a± = exp{iJ,/T) exp 



2k 



(B12) 



where ki = j/Tq. At zero temperature we used for the 
analytical continuation of the distribution function ||l8| 



fpik) = e{kF - \k\) 



(B13) 



with 



/^^(fc) = hm FA(fc) 

TT 





k + kp 




k — kp 




■| arctan 


A 


— arctan 


A 


} 



(B14) 



The calculation of Reej^e (B5) and Ime/je (B7) for com- 
plex frequencies will be d one analogously. The imaginary 
part of e/m in Eq. (B4) vanishes if A goes to zero and 
can therefore be dropped. 



APPENDIX C: THE THERMAL AVERAGING 

We define the thermal averaging by denoting that a 
function yl(e) should take its value onto the Fermi level 
for T = 

j de6{e-ef)A{e) 
1 



T 



defil-f)Aie) 



no J (2'iTy 



(CI) 



with the density of state nu — mpfT/ (27r^) for low tem- 
peratures. 



APPENDIX D: CALCULATION OF COLLISION 
INTEGRAL 



-A 



2^2 



rp6 



(D2) 



(D3) 



We perform the Fermi integral in (JS^) as above 
and have for Im 

00 

/m(±f^) ^ J "^^1 "^^2 dx3 dx4 S{Ax ± fl) 

-A 

where fJ = uj/T. Writing ( p^ ) in the form 

oc 

/m(±f^) ~ J "^^2 dx3 dx4 6{Ax ± n) 



-A 



fO fO 70 70 I f 2 .0 70 70 ( Tit _ 1 \ 

(D5) 



one gets for t he r esult (D2). Applying the stan- 
dard identity Q in (D5) we have 



rp3 r 

/„(±17) = — dxi 

-A 



(D6) 



so that the final result for l/^ in (D3 ) reads 
"2 



1 /W\2 



rf 



3 2 VT 
4 VttT/ 



(D7) 



Performing the collision integral about Fermi functions 
in (^) with the dimcnsionless variables x = {e — fj,) /T 
and A = /i/T we get 

00 

4 = T^j dx, dx, dx, dx4 5{Ax)flfmX^ (Dl) 



'here: A = f > 
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1. Without thermal averaging 



Calculating the Fermi integral in (^0|) wit hout the ther- 
mal averaging (app. |^) we start from (D4) and get 



/m(±i^) = r / dx2 dx3 dxiS{Ax ± fl) 



-A 



(fOfOfQ, rO 70 fO \ 
\J X2'' X3J X4 ~ J X2J X^J X^J 



(D8) 



Rewriting (D8) in the form 



(D9) 



/„i(±r2) — J dx2 dx3 dxi 5{/\x ± f2) 

-A 

we apply the standard identity [3^ ] and find for 



(DIO) 



Since e/ « ^ we have xi ^ and get for 7/^ ( |D3[ ) 
T2^2 r 



1 



The resulting relaxation time reads finally 



1 


3 1 


1 + 




Tm(w) 













(Dll) 



(D12) 
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